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ABSTRACT 

It is widely accepted that supersonic, magnetised turbulence plays a fundamental role for star 
formation in molecular clouds. It produces the initial dense gas seeds out of which new stars 
can form. However, the exact relation between gas compression, turbulent Mach number, and 
magnetic field strength is still poorly understood. Here, we introduce and test an analytical 
prediction for the relation between the density variance and the root-mean-square Mach num- 
ber ^ in supersonic, isothermal, magnetised turbulent flows. We approximate the density and 
velocity structure of the interstellar medium as a superposition of shock waves. We obtain the 
density contrast considering the momentum equation for a single magnetised shock and ex- 
trapolate this result to the entire cloud. Depending on the field geometry, we then make three 
different assumptions based on observational and theoretical constraints: B independent of p, 
B oc p 1 ' 2 and B cc p. We test the analytically derived density variance-Mach number relation 
with numerical simulations, and find that for B cc p 1 / 2 , the variance in the logarithmic density 
contrast, <x 2 n . = ln[l + b 2 ^ 2 Bo/(J3{) + 1)], fits very well to simulated data with turbulent 
forcing parameter b - 0.4, when the gas is super- Alfvenic. However, this result breaks down 
when the turbulence becomes trans-Alfvenic or sub-Alfvenic, because in this regime the tur- 
bulence becomes highly anisotropic. Our density variance-Mach number relations simplify 
to the purely hydrodynamic relation as the ratio of thermal to magnetic pressure Bq — > co. 
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Detailed knowledge about the statistical characteristics of the den- 
sity structure is of pivotal importance for many fields in astronomy 
and astrophysics. Probability distribution functions (PDFs) of the 
density have been introduced as a simple and robust measure of 
the one-point statistics for many applications, ranging from cos- 
mology, w here the Press-Schechter formalism was primarily es- 
tablished l lPress & Schechten, 1 19741) . to star formation and theo- 
ries of the initial mass function or the core mass function (e.g. , 
Fleckl Il982l; IZinneckeJ, Il984l; IPadoan. Nordlund. & Jonesl 1 1997 : 
Klessen & Burkert, 2000; Li et al„ 2004; Hennebelle & Chabrier, 



20081 12009; Padoan & Nordlun d kOlilT 



the root-mean-square (rms) Mach number in supersonic turbulent 
flow is a key ingredient for analytical models of the star forma- 
tion rate JKrumholz & McKed , [200l ; I Padoan & Nordlund . 1201 ll) . 
and for the stellar initial mass function or the core mass function 
jPadoan & Nordlund I2OO2I ; iHennebelle & Chabrien. [2008, 2009). 
In this framework, supersonic turbulence plays a fundamental role 
in deter mining the den sity and velocity stati stics of the inter stellar 
medium felmegreen & Scal o, 2004; McKee & O strikeil 12007b and 
controls stellar birth (Mac Low & Klessen, 2004). Conversely, the 
importance of magnetic fields in the star formation process is still 
incon c lusive, despite decades of r esear ch (Mouschovias & Ciolek, 
1999; McKee & Ostriker, 2007; ICrutcher. Hakobian. & Troland , 



In the star formation context, the relation between the width of 
the density PDF - the density variance or standard deviation - and 
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2009; Crutcher et al., 2010; Bertr am et al.l 120121) . Hence, the ques 



tion of how magnetic fields affect the density variance-Mach num- 
ber rel ation is stil l not clearly answered, despite the empirical find- 
ings of Ostriker et al. ( 200 1 ) and the analytical ansatz provided by 
IPadoan & Nordlund ( 120 nTT 



For purely hydrodynamical, supersonic, isothermal, turbulent 
gas, the relation between the density variance and Mach number 
has been identified and widely studied in numerical simulations 
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(e.g., Padoan et al., 1997; Passot & Vazquez-Semadeni, 



Federrath et al., 2008a; Federrath, Klessen, & Schmidt. 



Federrath et al., 2010; Price, Federrath. & Brunt, 2011). This rela- 



tion is commonly assumed to be linear, 



p/po 



bJK, 



(1) 



where cr 1 , is the density variance (to emphasise the density fluc- 
tuations about the meanp , it makes sense to express the density in 
terms of the density contrast p/po), b is a proportionality constant 
of order unity as explained in more detail below, and ^£ is the rms 
Mach number. Usually, the density contrast is written in terms of 
its logarithm, s = ln(p/p ). 

Several authors have noted that the PDF of the logarithm 
of the density contrast s - produced by supersonic turbulent 
flow of isother mal gas - follows a pprox i mately a lognormal 
distribution (e.g. Vazquez-Semadeni, 1994; Padoan et al., 
Passot & Vazquez-Semadeni, 1998; Nordlund & Padoan, 
KlessenL 12000; lOstriker, Stone, & Gammiel 1200 ll; lli et all 12003 
Kritsuk et all 120071 ; [Federrath et al., 2008b; Lemaster & Stonej 

gom 



1997; 



1999; 



Schmidt et al., 



2009 



Glover ct al.. 



Federrath et al.. 



200: .. .. _ _ _ . . 

20ld ; iPadoan & Nordlundl KOI ll ; ICollins et all 1201 ll ; iPrice et al 
201 lh. 



Ps ds ■ 



V^Rrf 



exp 



(s - s y 



lo-l 



d.v, 



(2) 



where the mean sq is related to the density variance by so = -cr 2 J2, 
due to the constraint of mass c o nserva t ion. Besides t h e em pirical 
findings of Vazquez-Semadeni 1 1994D. IPadoan et alj i 19971) , and 
IPassot & Va zquez-Semadeni ( 1998), there is no clear explanation 
for the shape of the PDF. From a mathematical point of view, a log- 
normal distribution is the result of independent random perturba- 
tions driven in a stationary system (Pope & Ching, 1993) as a con- 
sequence of t he cen t ral limit theorem (Vazquez-Semadeni, 1994; 
Padoa n et all ll997J ; iNordlund & Padoanl 1 1999; Federrath et al., 



20101) . The physical interpretation is that density fluctuations 



present at a given location are produced by successive passages of 
shocks with amplitudes independent of the local density. For a log- 
normal distribution, the density variance - given by Equation (TJl - 
is equivalent to 



ln[l+fc 2 



(3) 



The parameter b in Equations 10 and I0 is related to the ki- 
netic energy injection mechanism - the forcing F, which drives the 
turbulence. Federrath et al. (2008b) found that b = 1 for purely 
compressive (curl-free) forcing, VxF = 0, while b = 1/3 for purely 
solenoidal (divergence -free) forcing, VF = 0. In a follow-up study, 
iFederrath et al.l (I2010Q showed that b increases smoothly from 1/3 
to 1, when the amount of compressive modes, -F CO mp/(-Fsoi+.Fcomp) is 
gradually increased from to 1. For the natural mixture of modes, 
fcomp/(fsoi + ^comp) = 1/3, which is also the mixture of forc- 
ing modes used in all our numerical experiments here, they found 
b a 0.4, so we will later use that value for comparing our analytic 
model with numerical simulations. 

When magnetic fields are included, the density variance is sig- 
nificantly lower than in th e unmagnetised case for simulations wit h 
Mach nu mbers J( > 10 tostriker et all 1200 it IPrice et all 1201 lh . 
Recently, Padoan & Nordlund (201 1) provided an analytical ansatz 
for the hydrodynamical density contrast in supersonic, turbulent 



flow, which in turn follows the approach of iDvson & Williams! 
(1980) for obtaining the density contrast for strong adiabatic 
shocks, but extended to the magnetic case. Their cr s -^£ relation 
was, however, not tested with numerical simulations. 

The density PDF may or may not deviate from a log-normal 
form when other processes - like heat exchange and gravitation 
- are included. For example, when a non-isothermal equation of 
state is considered, the PDF still closely follo ws a log-normal dis- 
tribut ion over a range of densities (see e.g., iGlover & Mac LowL 
120071) . However, depending on whether the equation of state 
is softer or harder than isothermal, it might acquire power- law 
tails either at high or low densities (Passot & Vazquez-Semadeni, 
1 1998 1; IScaloetallll998l;IWadal l200ll ; iLi. Klessen. & Mac Lowi , 
120031; iMcKee & Ostrikeil 120071) . The density PDF also devi- 
ates from log-normal when gravity is included. In this in- 
stance, the PDF exh ibits a power-law tail at high densities 
JKlessed. 1200(3; IFederrath et all [ 2008a; K ainulainen et all 120091 ; 
ICho & Kid , 1201 ll : iKritsuk, Norman. & Wagnel 1201 ll) . In addi- 
tion, turbulent intermittency also leads to d eviations from the log - 
normal PDF in the wings of the distribution ( Federrath et alll2010l) . 
Consequently, the accuracy of the measurement of the density vari- 
ance, using Equation (O, may be compromised depending on the 
importance of the different processes involved in real molecular 
clouds. 

Here, we present an analytical derivation for the a s -.M re- 
lation in supersonic turbulent isothermal gas includi ng magnetic 
fields. Our results are i n qual itative agreement with lOstriker et all 
(2001) and Price et all OOlll) . however, here we present quantita- 
tive predictions and tests. The present work is organised as follows: 
In ^2]we describe the analytical approach made for the <t s —j& rela- 
tion. In this section, we start with the study of the density contrast 
of a single shock confined into a cubic box, and then we extrapo- 
late it to the whole cloud in 92.11 In 92.2l we propose three a s -.M 
relations given by three different assumptions of the behaviour of 
magnetic fields with density. We test these predictions with numer- 
ical simulations in $3] and conclude in §4. 



2 ANALYTICAL DERIVATION 

Our basis for obtaining the density variance-Mach number rela- 
tionship involves determining how the density contrast changes 
with the Mach number. The density variance <t p/po and the density 
contrast are related by: 



p/po 



p_ 

Po 



ll dV, 



(4) 



where p is the local density, p is the mean density in the volume, 
and V is the volume of the cloud. The density contrast is a measure 
of the density fluctuations in the flow, and therefore it is useful for 
identifying the disturbances that originate from shock fronts and 
compressions. 

2.1 Density contrast in magnetohydrodynamics 

Supersonic turbulence in the interstellar medium generates a com- 
plex network of shock waves (or simply shocks). When the veloc- 
ity of the fluid exceeds that of sound, it leads to the formation of 
shocks that are one of the most i mportant distinctive effec ts of the 
compressibility of the fluid (e.g. lLandau & Lif shitzl 1 1 987f) ■ 

In order to study the density contrast in a molecular cloud, 
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we first consider the physics of the discontinuity formed by a sin- 
gle shock front. We then generalise t he results to the ensemble of 
shocks confined in a cloud. Following Leaueux ( 2005), we describe 
a shock by choosing two control surfaces, one on either side of the 
discontinuity, and parallel to each other. Let us choose the shock 
surfaces as the reference frame, such that the control surfaces are 
stationary with respect to the shock. We also define the "parallel" 
direction as the one parallel to the flow of gas through the shock 
(i.e., perpendicular to the shock front). From the well known equa- 
tions of fluid dynamics, it is then possible to derive equations that 
expresses the conservation of matter and momentum flux for a mag- 
netised inviscid, neutral fluid: 



V||,lPl = v l2 p 2 , 



and 



Pi 



Y\ 



~2~ 



■-P2 



<2 >' 

+ — - + 

72 



A±,2 

2~~ 



(5) 



(6) 



respectively. In these equations, the subscripts 1 and 2 indicate the 
pre-and post-shock conditions, respectively. The velocity of the gas 
into the shock is vy, while c s is the adiabatic sound speed, y is the 
ratio between the specific heats and Vax is the Alfven velocity, de- 
fined here as v A± = B ± /(4np) i/2 , where B ± is the magnetic field 
perpendicular to the flow direction. The post-shock density is de- 
scribed by p2 . 

We now make two important approximations. First, as we 
wish to focus on the role of magnetic fields in determining the den- 
sity variance, we assume that the gas is isothermal, deferring con- 
sideration of non-isothermal effects to future work. Our assump- 
tion of isothermality implies that c sl = c s2 = C s and y t = y 2 = 1. 
Second, as we are considering an entire molecular cloud, we ap- 
proximate it as an ensemble of shocks. We assume that we can 
express the average pre-shock velocity in terms of the rms velocity 
Vo - hereafter, the subscript "0" indicates the volume averages - as 
vL = Trv 2 ,, where the factor b depends on the nu mber of degrees 
of freedom available for the compressive modes (Federrath et al., 
2008b). We also assume that the typical pre-shock magnetic and 
thermal pressures are just those given by volume averages over the 
total volume, allowing us to write them in terms of the volume- 
averaged density p and the rms Alfenic velocity v A0 . Similarly, 
we assume that the typical pre-shock density is simply the volume- 
averaged density. Making these assumptions, and introducing the 
ratio of the thermal pressure to magnetic pressure 



p 

1 man 



(7) 



case where the gas flows in a direction parallel to the field lines, the 
field strength will be the same on either side of the shock despite 
the jump in density, and the field strength then will be independent 
of density, i.e., a = 0. In the other extreme case where the field is 
oriented at right-angles to the gas flow, the shock jump conditions 
for magnetic flux freezing imply that B oc p, i.e., a = 1 . Meanwhile, 
compression of an isotropic field along all three spatial directions 
gives B oc p 2/3 . However, for our "average shock", we expect be- 
haviour that lies somewhere between < a < 1. By looking at 
observations and existing simulations, we can get some guidance 
as to what this intermediate behaviour should be. 

Observationally, Crutcher ( 1999) presented a study of the 
magnetic field strength in molecular clouds measured with the Zee- 
man effect. He fitted the results with a power law B oc p a and found 
that a = 0.47 ± 0.08. ICrutcher. Heiles. & Trolandl J2003b provided 
additional support for this result. More recently, [Crutcher et alj 
(2010) have presented a detailed compilation of Zeeman data based 
on a much larger number of measurements. They find that at num- 
ber densities n < 300 citT 3 , the data is consistent with a field 
strength that is independent of density, while at higher densities 
they obtain B oc p° - 65 * 005 . 

From a theoretical point of view, IPadoan & Nordlundl ( 11999b 
noted that their B distributions close ly match the observa tional 
scaling given by ICrutcherl d 19991) and ICrutcher et all d2003l) . B oc 
p l/2 , for high B in their high Alfvenic Mach number regime. 
iKim. Balsara. & Mac Low! 120011 ) also study the relationship be- 
tween B and p, and find that a ^ 0.4, albeit with large sca tter, 
especially at low densities. Additionally. iBaneriee et alj 120091) re- 
port that the magnetic field strength appears to scale in their simu- 
lations as B oc p 1/2 for number densities 10 2 < n < 10 4 cnT 3 , al- 
though with sign ificant scatt er around this value. On the other hand, 
iHennebelle & Peraulti (2000) found that the magnetic field does not 
necessarily increase with the density. Asides from these reports, if 
the magnetic flux is not conserved, but increases due to turbulent 
dynamo amplification during compression, a can become larger 
than the v alues quoted above, depending on the Reynolds numbers 
of the gas dSchleicher et al.Ll20ldrSur et alll201oKFederrath et all 
120111) . Thus, even if the gas is compressed only parallel to the field 
lines, turbulent tangling of the field can lead to a > during com- 
pression. 

Given the different possible relations between the magnetic 
field strength and the density, we consider three cases to include in 
Equation {8}: the two extreme cases, where B is independent of the 
density, and where B oc p, and an intermediate case with B oc p [/2 . 
We also note that if we were to take instead the relation B oc p 065 
suggested by the most recent observational data, then we would 
obtain results quite similar to the B oc p 1/2 case. 



we can rewrite Equation (|6) as 



^ 2 H-^) + ^(i + /V) = (i + pV), 

P2 \ P2 / P2 V ' V ; 



(8) 



where the rms Mach number is given by ./£ = Vo/c„. 

In order to solve this equation for the characteristic density 
contrast associated with the shocked gas, p 2 /po, it is necessary to 
determine B 2 , the post-shock ratio of the thermal to magnetic pres- 
sures. The value of this will depend on the change in the magnetic 
field strength through the shock, which in turn depends on the ori- 
entation of the field with respect to the flow of gas through the 
shock. Using magnetic flux and mass conservation during com- 
pression, one can show that B oc p a with < a < 1, depending 
on the field geometry and direction of compression. In the extreme 



2.7.7 First case: B independent of p 

We start by considering one extreme, the case where B is inde- 
pendent of the density. In this scenario, Equation ([8) becomes a 
second-order equation, independent of the magnetic field strength 



\Po) v '\Po 



j + b 2 J( 2 = 0. 
This equation results in a density contrast 



P -± = b 2 ,, 
Po 



(9) 
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Equation ^ matche s the density con trast for the non- 
magnetic regime (see e.g. JPadoan et alll 19971) . This is not surpris- 
ing, because in this case we are assuming that the gas and the mag- 
netic field are not coupled. Therefore, amplification of the magnetic 
field with density is not expected under these conditions. 



2.1.2 Second case: B oc p l/2 

In the intermediate case in which B oc p 1 ' 2 , we again find a second- 
order equation for the density contrast, but with a dependence on 
the magnetic field expressed in terms of 8 . From Equation l[8}, we 
obtain 



(l +y3 ')(-) ~{b 2 ^ 2 + 1 +&')(-) +b 2 ^ 2 = 0. 
This equation has the solution: 



£ = b 2 Jt 2 
Po 



B + l 



(10) 



In other words, the effect of the magnetic field in this case is to 
reduce the density contrast by a factor Bq/(Bq + 1). We see from 
this that in the weak field limit where B — > <x>, we recover the 
hydrodynamical result, while for strong fields we have a smaller 
density contrast in the MHD case than in the non-magnetic case. 



2. 1.3 Third case: B oc p 

Finally, we investigate the other extreme case, where the magnetic 
field strength is proportional to the density. In this case, Equation 
([8} results in a third-order equation, 



« , (S) ,+ teJ-c +1+ « , )(S) +t "- a 

The solution for the density contrast is 

8o+^{\+B Q f + Ab 2 ^ 2 B \. (11) 



-4h 

Po 2 



2.2 Density variance-Mach number relation 

In the previous section, we presented three different expressions for 
the density contrast. They correspond to three different assumptions 
regarding the relationship B oc p", with a = 0, 1/2, and 1. We now 
determine the density variance of a fluid in which there are many 
shocks, for each of these three cases. 

We start by noting that in a highly supersonic flow, the dom- 
inant contribution to the integral in Equation (4) will come from 
shocked regions, and thus we can consider this equation as a vol- 
ume average over an ensemble of many shocks. We next assume 
that we can approximate the value of this integral with the result 
of integrating over a single "average" shock of the kind considered 
in the previous section. As we already know the density contrast 
of this representative shock, the only thing that remains to be done 
before we can solve Equation (4) is to determine the appropriate 
volume over which to integrate. 

We approximate the cloud as a cubic box of side L, and con- 
sider an infinitesimal part of its volume dV that encloses one shock. 
Therefore, the size of dV depends on the size of the shock itself 



dV ~ dV, 



si,- 



(12) 



To define the shock volume, we make us e of an approxima- 
tion introduced by IPadoan & Nordlundl ( 1201 II) , where the volume 
of the shock is given by the area of the box face times the shock 
width A, V s h = L 2 A. However, in the absence of viscosity, it is not 
straightforward to define th e shock width A. Therefore, we follow 
IPadoan & N ordlund ( 20 1 1 ) and assume that the shock width, if the 
compression is driven at the box scale, is given by 



A - 6Lp /p 2 , 



(13) 



where 6 is the integral scale of the turbulence. Then, the volume of 
the shock V^ is given by 



Pi 



(14) 



For turbulence driven on large scales, as appears to be the case in 
real molecular clouds l lOssenkopf & Mac Lowi |2002| ; iBrunt et all 
2009), we have 0^1. Having made the assumption that the ap- 
propriate volume over which to average is the volume of our repre- 
sentative shock, and considering Equation [15] we approximate dV 
by 



w=i? md p2 



\P2) \Poj 

Finally, inserting Equation (Q3) into Equation ©, yields 



(15) 



plpa 



!:( 



-2 In 



(16) 



P0 d |P2 _ P_ P0 

Pi) \po) Po P \P0) 

It is important to note that in this formulation, Equation d!6t is 
physically meaningless if the lower limit of the integral is set be- 
tween < p/po < 1. It is due to the definition adopted for the 
shock width (Eq.ll4ll. where the shock thickness is defined only for 
P2/P0 > 1- For highly supersonic turbulence, which is the regime 
that concerns us, the assumption p » po is valid. Then, the first 
term in Equation dl6) dominates the variance and we get 



p/po 



p_ 

Po' 



(17) 



For practical reasons, we prefer to consider the variance of 
the logarithm of the density contrast, s = ln(p/p ), instead of the 
variance of the linear density when we will compare this analytical 
model with numerical simulations. These varia nces are related by 
(e.g jFederrafh et alll2008bl;|Price et alll201 lh 



ln[l 



'p/poj 



(18) 



We now insert the three cases considered in 92. II into Equa- 
tion dl8K in order to obtain the density variance-Mach number re- 
lation. The subscripts of the following results are chosen based on 
the value a = 0, 1/2 and 1 of the B oc p a relationship. 

• B independent of p 

The density variance in this case is exactly the same as for the 
purely hydrodynamical, isothermal case, 

cr 2 = ln[l+/r.^ 2 ]. (19) 
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• B oc p 1 ' 2 

In this case, the density variance is: 



s,l/2 



In 



A, 



(20) 



\Po + 1 , 

This relation is similar to Equation d!9t except for a correction 
factor due to the effects of magnetic fields, which is a function of 
the plasma /3 only. 

• B oc p 

Finally, the density variance-Mach number relation in this case 
is given by 



In 



i + \ f-i -A +Vo+A0 7 +4R 



^ 2 A, 



(21) 



The density variance has a strong dependence on /3o, leaving the 
rms Mach number as a marginal quantity in this relation. 

In the last two cases, when p — > 0, the Alfenic velocity is 
much higher than the sound speed, and both relations approach 
zero. In this scenario, the magnetic pressure is infinitely large and 
prevents density fluctuations from forming. The gas is "frozen" in 
the magnetic field. In the opposite limit, wheny8 — > °°, Equation 
d20t and Equation d21t simplify to the purely hydrodynamical case, 
as expected. In the next section, we are going to test these cases 
with numerical simulations. 



3 NUMERICAL TEST OF THE ANALYTICAL MODEL 

3.1 Simulations 

We have performed simulations of the evolution of the turbu- 
lent, dense, inviscid, magnetised (MHD) and unmagnetised (HD), 
isothermal interstellar mediu m using a mod i fied version of th e 
zeus-mp hydrodynamical code (Norman, 2000; Haves et al., 2006). 
We neglect chemical reactions in order to study the effects of mag- 
netic fields in molecul ar clouds, l eaving the inclusion of the effects 
of chemistry (Gl over et all 120101) for a future study. 

Each of our simulations begins with an initially uniform gas 
distribution, with a mean hydrogen number density of n = 1000 
citT 3 and a resolution of 256 3 cells. The initial velocity field is 
turbulent, with power concentrated on large scales, between wave 
numbers k = 1 and 2 and with an initial rms velocity of 5 kms" 1 . 
Moreover, we drive the turbulence so as to maintain approximately 
the same rms velocit y throughout the sim ulati ons, following th e 
method described in lMac Low et alj dl998h and lMacLowl J 19991) . 
We do not perform a Helmholtz decomposition of the force field, 
and thus the turbulent forcing consists of a natural mixture of 
solenoidal and compressive modes , i.e., F sol /(F 50l + F comp ) » 2/3. 
Note that iFederrafh et al.1 J2008bl l2010h tested the two limiting 
cases of purely solenoidal (divergence-free) and purely compres- 
sive (curl-free) forcing, as well as various mixtures of solenoidal 
and compressive modes of the turbulent forcing. They found a 
strong influence on the density PDF, producing a three times larger 
standard deviation for compressive forcing compared to solenoidal 
forcing. They parameterised the influence of the forcing by intro- 
ducing the ^-parameter in Equation ((3}. Purely solenoidal forcing 
is characterised by b = 1/3, while purely compressive forcing gives 
b = 1, For the natural mixture, they find b w 0.4. Using the present 
set of numerical models, we confirm that using b = 0.4 for the 
natural mixture of forcing modes used here gives the best fits with 



our analytically derived density variance-Mach number relation. 
The temperature of the gas is constant and fixed to an initial value 
T = 1062, 170, 42 and 15 K, in order to sample a large set of 
Mach numbers {^£} ^ 2, 5, 10 and 17, respectively. We adopt 
periodic boundary conditions for the gas using a cubical simula- 
tion volume with a side length L = 20 pc, such that the turbu- 
lent crossing time, r cros s = L/(2c s ^) x 2 Myr. We present re- 
sults from t = 3 r cross » 5.7 Myr, sampled every 0.17 7" cross , and 
evolved until t = 4r cross as 7.6 Myr. This period of time is long 
enou gh to expect the turbulence to have reached a statistical stead 



I 



state (JF ederrath, Klessen, &JchmidtL[2009; Fed errathet^ll2oT 
Glover et al., 2010; Price & Federrath, 2010). This simulation time 
might be also short enough to obtain reliable results for the initial 
phase of star formation, when self-gravity did not yet have a large 
effect on the dynamics. In order to concentrate on turbulent com- 
pression alone, we neglect self-gravity in the present experiments. 

For the MHD cases, the simulations begin with a uniform 
magnetic field that is initially oriented parallel to the z-axis of the 
simulation. Four of these simulations begin with an initial magnetic 
field strength 6, = 5.85 /iG, which is our standard magnetic field 
strength hereafter. We also perform three MHD runs with B, = 10, 
20 and 60 /jG, with ^ = 10, to check the behaviour of the re- 
sults with increasing magnetic field strengths. We note that as the 
simulations run, dynamo amplification can lead to increased field 
strength, and thus we use the instantaneous magnetic field strength 
to compute y8 . Nevertheless, for simplicity we use the initial value 
of the magnetic field strength to label runs MHD-B2, MHD-B20 
and MHD-B60. 

In Table 1 , we list the simulations that we have performed. In 
our labels, we use "H" to denote a hydrodynamic run and "MHD" 
to denote a magnetohydrodynamic run. Our multiple runs with 
fixed (or zero) magnetic field strength but different sound-speeds 
are labelled with an "M", followed by the (approximate) rms Mach 
number of the simulation. Finally, the three runs in which we ex- 
amined the effect of varying the initial magnetic field strength are 
labelled with a "B", followed by the initial field strength in fiG. 
In Table 1, we also list the values of the quantities: (5q, the rms 
Alfvenic Mach number ./id.o = vo/v^.o and the sonic Mach num- 
ber. They are measured in every cell and then are spatially averaged 
over the datacube. The brackets denote the time average over the 
seven snapshots, and the lcr shows the temporal standard deviation 
around the mean values. 



3.2 Statistical Analysis 

In this subsection, we explain the method used to measure the den- 
sity variance for every snapshot in our si mulations using the PDF as 
a robust statistical tool for this analysis l lPrice et all 120111) . Then, 
we parameterise the instantaneous /?o in terms of ./£ ', in the direc- 
tion of testing numerically the a-./£ relations presented in 32.21 
Finally, we present the comparison between our analytical model 
and the simulations. 



3.2.1 Probability Density Function (PDF) 

In Figure [T] we plot the volume-weighted dimensionless density 
PDFs for MHD and HD isothermal gas with the same Mach num- 
ber for comparison. For these simulations, we find that all the PDFs 
have a log-normal shape around their peak. However, the PDFs de- 
viate from log-normality especially in the HD simulations at low 
densities, being more evident for ./£ > 5. The error bars in this 
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Table 1, List of simulations. 



B, 


(Pa) ± lo" 


<^A,o) ± lo" 


<o-j> ± io- <./#> ± lo- 


HD-M2 


CO 





om ± 0.02 2.21 ± 0.02 


HD-M5 


oo 





1.3 ±0.1 5.4 ±0.1 


HD-M10 


oo 





1.7 ±0.1 10.6 ±0.2 


HD-M17 


oo 





1.92 ±0.09 17.6 ±0.5 


MHD-M2 5.85 


25 ± 5 


8.1 ±0.9 


0.69 ± 0.02 2.09 ± 0.02 


MHD-M5 5.85 


4.8 ± 0.4 


8.4 ± 0.8 


1.18 ±0.04 4.98 ±0.07 


MHD-M10 5.85 


1.4 ± 0.5 


9±3 


1.47 ±0.06 10.2 ±0.3 


MHD-M17 5.85 


0.3 ± 0.1 


7±2 


1.61 ±0.06 16.8 ±0.5 


MHD-B2 2 


11.3 ± 0.5 


27 ±2 


1.58 ±0.09 10.5 ±0.2 


MHD-B20 20 


0.083 ± 0.005 


1.94 ±0.06 1.48 ±0.01 9.9 ± 0.2 


MHD-B60 60 


0.030 ± 0.001 


1.24 ± 0.03 1.34 ± 0.01 10.3 ± 0.1 



Bj - initial magnetic field strength in /iG. 

Bo - mean thermal to instantaneous magnetic pressure ratio. 

./sCi.o - rms Alfvenic Mach number. 

<t s - density variance. 

./# - rms Mach number. 

The brackets indicate the time average calculated over the snapshots after 

averaging over the spatial coordinates. 

figure show the lcr variations around the time average. We see 
that these variations cannot explain the tail at low densities. There- 
fore, this deviation is not explained by intermittency fluctuations, 
and deserves further study. However, the low-density tail does not 
significantly affect our a s estimates, because the variance is com- 
puted from a log-normal fit in a limited inte rval around the pe ak, 
giving the most reliable estimates of cr s (see lPrice et all 1201 If) . In 
this sense, the trend of the time averages observed between MHD 
and HD simulations shows the magnetic field acting as a density 
cushion, preventing the gas from reaching very low densities dur- 
ing local expansion. As a consequence, there are larger parts of the 
volume with density p « po in the MHD case than in the HD case. 
In order to avoid contamination from intermittency, numerical 
artefacts, etc., in the wings of the PDFs, we perform a Gaussian 
fitting only in a data subset selected by s, in each simulation. This 
subset consists of 60% of the number of bins considered to calcu- 
late the density PDF which are distributed symmetrically around 
the mean, sq. Then, we fit the Gaussian profile given by Equation 
10 to obtain cr s in every snapshot of the simulations. 



3.2.2 Density variance— rms Mach number test 

In the interest of comparing the density variance-Mach number re- 
lation, given by Equation J20b and Equation \2\\ , with the results 
obtained in the previous subsection, we parameterise the thermal- 
to-magnetic pressure ratio in terms of the rms Mach number for our 
sequence of simulations. In this sense, we rewrite Equation 10 as 



Po 



',1,0 



(22) 



Note that this parameter is calculated considering the instantaneous 
magnetic field strength and not the initial value. 

Next, we select the four MHD simulations with different rms 
Mach number, but the same initial magnetic field strength, and 
use a linear regression considering the logarithm of Equation {22): 
log 10 /3o = log 10 C — 2 log 10 ./£ . From the fit shown in Figure|2] we 
find C = 1 1 1 ± 4. In Figure [2] we plot Bq as a function of the rms 
Mach number for the different snapshots. The triangles show [} a for 






1 


10" 


-1 


10' 


■2 


io- 


■3 



-10 



-5 

s=ln(p/p ) 



Figure 1. Dimensionless density PDF for magnetised and unmagnetised 
molecular clouds with the same initial conditions, no = 1000 cirT 3 , and 
same turbulent rms velocity, but different sound speed. The most significant 
features are: 1) the density variance increases with Mach number, and 2) the 
density variance decreases with magnetic field strength. These simulations 
have a ratio between thermal pressure and magnetic pressure Bq < 10. All 
simulations have a resolution of 256 3 zones. 



the selected simulations with (,/#> » 2, 5, 10 and 17, while the 
curve shows the linear regression. 

In Figure [3] we combine the dimensionless standard deviation 
cr s , obtained from the fit over the numerical PDFs for every snap- 
shot, and the analytical prediction for the three cases of B oc p a 
- with a =0,1/2, and 1 - as a function of the rms Mach number. 
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Figure 2. Parameterisation of /3o = ^th/^o.raaa \vith respect to the rms Mach 
number for the subset of simulations with roughly constant Alfvenic Mach 
number, ./^a.o ~ 8 (see Table [TJ. The curve is a linear regression of the 
MHD simulations with B, = 5.85/iG. The linear regression performed to 
the logarithm of Equation l(22j gives /J = (1 1 1 ± 4)„#~ 2 . 



For the triangles around a given {./£}, the HD simulations exhibit 
larger cr s compared with the MHD simulations, as was expected 
from Figure [TJ For comparison, we plot the analytical prediction 
given by Equation J 19b. ov.q. T his result matches the prediction 
provided bv lPadoanet al.l J1997T) . However, instead of using their 
roportionality parame ter b « 0.5, we used the input value b = 0.4 
Federrath et al.l 120 lOl dashed line), which is the result of the nat- 
ural mixing of solenoidal and compressive modes in the turbulent 
forcing field. We also plot the two extreme cases for the unmag- 
netised gas, <t sHD , with b - 1/3 (lower dotted line) for purely 
solenoidal forcing and b = 1 for purely compressive forcing (upper 
dotted line) for comparison. 

In the same Figure, we superpose Equation J20I light grey 
solid line) and Equation J21I dark grey solid line), both again with 
b = 0.4. We find than the best agreement with the MHD simula- 
tions is given by Equation J20b . that is cr,,i/2- The result obtained 
for the first case - B independent of density (Equation 1 19ll - may 
account only for low Mach number zones. This c ase might be ap- 
propriate for diffuse clouds iCrutcher et al.luOlOO , where the mean 
sound speed of the cloud may be of the same order as the rms ve- 
locity. Here, at ^£ ~ 1, all the three cases converge to the HD 

result. 

Our re sults are qualitativ ely in agreement with lOstriker et alj 

feOOll) and |Priceet~al] < l201ll) . These authors find that the density 
variance in magnetised gas is significantly lower than in the HD 
counterp arts for simulations wi th a Mach number ..# > 10. In 
addition, ICho & Lazarianl I I20031) study the density contrast result- 
ing from the Alfvenic waves, slow and fast magneto-sonic waves 
originating in different environments. The authors concluded that 
the three kinds of waves can coexist in those environments. In the 
regime that concerns us, B » 1 and 5 < ./# < 10, their density 
contrasts closely match ours. 

To test the validity of our results for different Alfvenic 
Mach numbers, we also performed three simulations with an ini- 
tial magnetic field strength different from the standard one, with 
^A,a ~ 27, 1.9, and 1.2, at (^£) » 10 (empty squares in 
Figure [3}- Our model works well for ^Ci_o <: 6, but breaks 
down for our test with X < 2. The break occurs when the 



turbulence becomes trans-Alfvenic or sub-Alfvenic, i.e., when 
Sa,o 5: 2. This is due to anisotropics arising in this case, i.e., 
the turbulence is no longer isotropic, as can be seen in Figure 
[4] This is because the back reaction of the magnetic field onto 
the flow is extremely strong for flows perpendicular to the mag- 
netic field lines, if the turbulence is trans-Alfvenic or sub-A l fvenic 
(see e.g.,|C ho & Lazarian, 2003; IBrunt. Federrath. & Price! 120101 : 
lEsquivel & Lazarianl 1201 11) . Since our analytic derivation is based 
on an ensemble average (Eq. |4). assuming statistical isotropy, the 
anisotropies are the most likely cause for the limitation of our 
model to super- Alfvenic turbulence. In Figure [5] we show our pre- 
diction (Eq. I20I L I for a fixed Mach number ..# x 10 and forcing 
parameter b » 0.4, which fits very well the data with ,/£a,o <; 6. 
These simulations show high dispersion - around the time aver- 
age - in the density variance and the rms Alfvenic Mach number 
showing the fluctuations of the gas caused by the turbulence dom- 
inating the dynamics of the flow, in contraposition of the simula- 
tions with small Alfvenic Mach number. In the same Figure, we 
also plot the model curve Eq. d20t for the same sonic Mach number 
10 and b = \. Although our turbulent forcing in t he simulations is 
by definition mixed, and thus we expect b w 0.4 JFederrath et all 
2010), we find it interesting to note that b = 1 - corresponding 
to purely compressive forcing - gives a good fit to the data with 
very low Alfvenic Mach number, ^#a,o ^ 2. We speculate that the 
density field for very high magnetic field strengths and thus very 
low Alfvenic Mach number starts behaving as if it was driven by 
purely compressive forcing. This is very different from the com- 
pression obtained with solenoidal or mixed forcing, but more simi- 
lar to compressiv e forcing, which also directly compresses the gas 
JFederrath et al.l ;2008b). More data at ^a,o :£ 2 would be needed 
to sample this region and the transition from b = 0.4 to 1 in detail, 
and we just note here that b = 1 seems to provide a good fit for 
^a,o :£ 2, given the data at hand. 



4 CONCLUSIONS 

We presented an analytical prediction for the density variance- 
Mach number relation in magnetised supersonic turbulent gas. In 
this formulation, we considered three different cases for the relation 
between the magnetic field strength and density. The first case as- 
sumes that B is independent of p, the second assumes that B oc p 1/2 , 
while the third is given by B oc p. The three resulting cr s -./£ rela- 
tions were tested against numerical simulations. From this analysis 
we conclude that: 

• If 6 is independent of the density, we recover the hydrody- 
namical prediction of Padoan et al. ( 1997). In this case, the gas and 
the magnetic field are not coupled. Therefore, an amplification of 
the magnetic fie ld with the shock is not expected. Observationally, 
ICrutchen (1 19991) found that the magnetic field was independent of 
the density for diffuse clouds, corresponding to low rms Mach num- 
bers, ./# < 2. In this regime, all our predictions converge to the 
purely hydrodynamical a s ~M relation. 

• For the second case, B oc p 1/2 , we found a one-to-one re- 
lation between ^#, /3o and the density variance. This cr,-^# re- 
lation (Eq. |20t matches very well our numerical test considering 



1 Equation J2Q) has been written in terms of the instantaneous Alfvenic 
Mach number (Eq. I22t . yielding the relation for the density variance: 



o-; i/2 =inU + 2b 2 y/ 



K\.u 



.4.0 " 



')]. 
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Figure 3. Standard deviation of the dimensionl ess density contrast, plotted as a function of the rms Mach number. Circles show the purely hydrodynamical 
simulations that follow very well the IPadoan et al.l 419971) prediction, u sHD = ln(l + b «# ), with b = 0.4, expected for mixed-mode turbulent forcing 
iFeder rath et alll201(j dashed line). The d otted lines are for comp arison with purely hydrodynamical model, assuming b = 1/3 for purely solenoidal forcing 
and b = 1 for purely compressive forcing (Federrath et al., 2008b). Triangles show the MHD simulations and the two formulas, Eqs. )20t and 12 IK obtained 
in this work: o" s- i/2 = {ln[l + b 2 ^ 2 fSo/(fio + 1)]} (light g rev solid line), and <x, i (dark grey solid line). Those curves are plotted for b = 0.4, and using 
our parameterisation, y8o = (111 ± 4)^~ 2 from Fig. [2] Squares, stars and diamonds show the additional MHD simulations with different rms Alfvenic Mach 
number, Jf. Afl as 27 (B, = 2/iG), ^ A ,o ~ 1.9 (B, = 20 //G), and J( Afs x 1.2 (B, = 60pG). 



(\=11.2 tJ r=10.a //,=25 |3 = 1.0. /^9.9, /4=7.3 p o = 0.09 t J/ =9.7 ff K =2 f3 o =O.03 t //=10.5 %1.7 




Figure 4. Density slices of the simulations at t = 6 Myr. The mean magnetic field is oriented along the vertical axis. From left to right: initial magnetic field 
strength B, = 2, 5.85, 20 and 60 pG. The turbulence remains isotropic for super- Alfvenic gas ^a,o * 1. but when it becomes trans-Alfvenic or sub Alfvenic 
(^A,o £ 3), the turbulence becomes highly anisotropic. 
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Figure 5. Standard deviation of the dimensionless density contrast, plotted 
as a function of the instantaneous rms Alfvenic Mach number at (^#) a 
10. The different symbols show snapshots of simulations with ./#a,o time 
averages: (^A,a) ~ 27 (squares), <.-#a,o> ~ 9 (triangles), (./#a,o) ~ 1-9 
(stars), and <^#a,o) ~ 1-2 (diamonds). When the turbulence becomes trans- 
Alfvenic or sub- Alfvenic, (.$a,q) 5> 2 (stars and diamonds), anisotropics 
arise in the gas, because the back reaction of the magnetic field onto the 
flow is extremely strong for flows perpendicular to the magnetic field lines. 
The grey curve shows our prediction ' using b a 0.4 that fits very well the 
data. Meanwhile, the black curve shows our prediction 1 considering b = 
1 (corresponding to purely compressive forcing). Although our turbulent 
forcing in the simulatio ns is by definition mixed, and thus we expect b a 0.4 
IFederrath et aUfeOlOl) . it is noteworthy to say that b = 1 gives a good fit to 
the data with very low (,J&a,o) :£ 2. 



b = 0.4, which is the input for the natural mixture of compressive- 
to-solenoidal modes in the turbulent fo rcing field. This resu lt is 
in agreement with the ones presented bv lOstriker et al.l (120011) and 
IPrice et al.l ( 1201 ID . where they found lower <j s than in the un mag- 
netised case for ^# > 10. Moreover, Icho & Lazarianl (120031) pre- 
sented a density contrast that closely matches our result for /?o ~ 1 
and 5 < J( < 10. 

• For the last case, B oc p, the a s -./M relation (Eq.l21t predicts 
a lower density variance than measured in our numerical simula- 
tions for ./# > 5, because our simulations are closer to B oc p 1/2 . 
However, the relation given by Equation l !21t would fit better, if 
B oc p. 

• The a s -^/£ relation obtained for B oc p [/2 works very well for 
intermediate to high Alfvenic Mach number, ^a,o ^ 6, but breaks 
down for ^a,o S= 2 at (^) » 10. This probably occurs because in 
the presence of strong magnetic fields, the turbulence is no longer 
isotropic. This is because the back reaction of the magnetic field 
onto the flow is very strong for flows perpendicular to the magnetic 
field lines. 

Magnetic fields act as a density cushion in turbulent gas, pre- 
venting the gas from reaching very low densities as well as very 
high densities. We conclude that magnetic fields are an important 
mechanism for shaping the density variance-Mach number rela- 
tion, and therefore will change the quantitative predictions in mod- 
els of the star formation rate, initial mass f unction or core mass 
function that depend on these quantities (e.g. Krumholz & McKee, 
2005 1; IPadoan & Nordlundl l201lL |2002| ; iHennebelle & Chabriei , 
2008J I2009). 
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